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2.6 EL PRINCIPIO DE LOS TRABAJOS VIRTUALES EN EL MÉTODO DE 
EQUILIBRIO 


En estas lecciones se ha utilizado previamente el PRINCIPIO DE LOS TRABAJOS VIRTUALES (P.T.V.) en 
una de sus aplicaciones conocida por el alumno, el TEOREMA DE LA CARGA UNIDAD. En este teorema, 
como movimientos y deformaciones, compatibles entre sí, se toman las del “estado real” de la estruc- 
tura, y como “estado virtual” se toma un sistema de acciones y solicitaciones, equilibradas entre sí, 
elegido adecuadamente para el propósito de calcular un determinado movimiento del “estado real” de la 
estructura. 


Posteriormente, en el apartado 2.4 se han formulado, como ejercicio puramente algebraico, los trabajos 
virtuales externo (We) e interno (Win) resultantes de adoptar como “sistema estático” las acciones y 
solicitaciones (f, s) del “estado real” de la estructura, que naturalmente deben cumplir las ecuaciones 
de equilibrio del problema (£=H-s); y como “sistema cinemático” los movimientos y deformaciones (u, 
e) de un “estado virtual” elegido arbitrariamente, con la única condición de cumplir las ecuaciones de 
compatibilidad (e=B-u). Se comprobó la igualdad Wer: =W ¡n; en el caso particular considerado. 


Allí se introdujeron las expresiones W ey =f-u y Winy=s'-e, que por ser f, s, u, e “vectores” (matrices 
de una sola columna), simplemente expresan que el trabajo virtual externo (We,:) es la suma Y f;-u; 
(suma de los productos de cada término del vector de acciones por el correspondiente término del vector 
de movimientos), y que el trabajo virtual interno (Wins) es la suma > s;-ej (suma de los productos de cada 
término del vector de solicitaciones por el correspondiente término del vector de deformaciones). Estas 
expresiones son válidas, lógicamente, sólo si los términos f; y u; se ordenan de manera “correspondiente” 
(desde el punto de vista del trabajo mecánico) en f y u; y análogamente para los términos de s y e. 


Puesto que las expresiones W..,¿=f'-u y Win;=s'-e son básicas para la validez de las demostraciones 
incluidas al final del presente apartado, se deduce la regla de “correspondencia” en la ordenación de los 
términos de los vectores antedichos, y por lo tanto de las ecuaciones que los ligan, regla que se introdujo 
en un aparado precedente y que ahora se justifica. 


La expresión Went=f!-u= Y £;-u, permite generalizar el concepto de “acción” (f;) a cualquier magnitud estática, o 
conjunto de ellas, que “trabajan por” efecto de un cierto movimiento u;. Esta generalización es importante en el análisis 
de pórticos desplazables de edificación, para los cuales, tal como se verá en apartados posteriores, frecuentemente se ligan 
determinados movimientos (por ejemplo, el desplazamiento horizontal u¿ de todos los nudos de una misma planta), y como 
consecuencia de ello la acción f; puede incluir varias fuerzas aplicadas a la estructura. 


Por otra parte, en los apartados previos dedicados a estructuras intraslacionales, se ha encontrado, por 
formulación directa de las matrices de equilibrio (H) y de compatibilidad (B), que ambas son transpuestas 
entre sí (H=B*; B=H?). 


En el presente apartado se demuestra, aplicando el PRINCIPIO DE LOS TRABAJOS VIRTUALES, que dicha 
propiedad es general (dentro de la teoría de primer orden), lo que tiene una enorme trascendencia en la 
aplicación práctica del método de equilibrio. En efecto, como consecuencia de dicha propiedad (H=B*; 
B=2H?), no es necesaria la doble formulación explícita de las ecuaciones de equilibrio (f=H-s) y de las 
ecuaciones de compatibilidad (e=B-u) de cada problema, puesto que de uno de los sistemas de ecuaciones 
(por ejemplo, el de las ecuaciones de equilibrio), se infiere el otro (las ecuaciones de compatibilidad; o 
viceversa). 


En las presentes lecciones se considera como objetivo docente prioritario saber formular explícitamente 
las condiciones de equilibrio de cada problema, como forma didáctica de continuar, incluso para las es- 
tructuras hiperestáticas, trabajando sobre las condiciones de equilibrio en tanto que herramientas funda- 
mentales del análisis de estructuras, cuya estabilidad depende ante todo de que no se vulnere ninguna de 
dichas condiciones. 


La relación B=H! es una alternativa obvia para hacer de las ecuaciones de equilibrio las fundamentales 
del análisis de estructuras hiperestáticas de barras, ya que éstas pueden formularse con sencillez e incidi- 
endo en conceptos que se han puesto en práctica previamente en el análisis de las estructuras isostáticas. 
Una vez escritas las ecuaciones de equilibrio, de la matriz H correspondiente se infiere automáticamente 
la matriz de compatibilidad B por una simple transposición de filas por columnas. 
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De esta manera no es necesario, ni tan siquiera, formular explícitamente las ecuaciones de compatibil- 
idad e=B-+u; pero a lo largo de los diversos ejemplos se insiste en la necesidad de saber “interpretarlas” 
(se escriben las diversas filas de e=B-u y se verifica su coherencia geométrica sobre un dibujo de la 
deformada de la estructura), como procedimiento eficaz de: 


e Comprobar que la matriz B, y por lo tanto también la H, son correctas. 


e Comprender cómo se mueve y deforma cada tipo estructural. 


Lógicamente es perfectamente válido partir de las ecuaciones de compatibilidad e=B-u, extraer de 
las mismas la matriz B, y obtener H por transposición (o bien, directamente formar la matriz de rigidez 
K=B*DB), siempre que se compruebe, dada su absoluta importancia, que las expresiones f=B*.s son 
las ecuaciones de equilibrio correctas del problema. 


La simetría de las ecuaciones de equilibrio y de compatibilidad otorga a las mismas un carácter absolutamente dual en 
el formalismo de los métodos de equilibrio, aunque por razones diversas puede reforzarse más o menos el estudio de una 
de ambas familias de ecuaciones. Estas lecciones, aprovechando la formación previa del alumno en estática, y primando la 
enseñanza del equilibrio como condición fundamental de la seguridad estructural, se decantan por tomar las ecuaciones de 
equilibrio como ecuaciones de partida para el análisis hiperestático. 


El estudiante de cursos de especialidad encontrará que en el cálculo de medios continuos por el método de elementos 
finitos, el método de equilibrio se aplica a la inversa: se formulan las ecuaciones de compatibilidad (si bien de manera sólo 
incompleta por utilizarse campos aproximados de movimientos, con un número finito de parámetros); y a partir de las 
mismas se infieren las de equilibrio (que consecuentemente tienen carácter aproximado). La expresión típica de la matriz 
rigidez de cada elemento finito responde al formato B*DB, con algunas variantes como la de una integración extendida al 
área O volumen del elemento. 


Como consecuencia de la relación H=B*, dentro del análisis elástico lineal, en que la matriz consti- 
tutiva es siempre simétrica (D=)D*), también lo es la matriz de rigidez K=B*DB, como se demuestra 
posteriormente. Las implicaciones prácticas de esta propiedad son importantes, ya que permite simpli- 
ficaciones en el almacenamiento en memoria de dicha matriz y en los algoritmos de resolución de las 
ecuaciones Ku=f 


A continuación se recuerdan las condiciones de aplicación del PRINCIPIO DE LOS TRABAJOS VIRTUALES: 


e Las solicitaciones y las acciones deben estar equilibradas (SISTEMA EQUILIBRADO) 
e Las deformaciones y los movimientos deben ser compatibles (SISTEMA COMPATIBLE) 


e El SISTEMA EQUILIBRADO y el SISTEMA COMPATIBLE pueden no tener relación de causa a efecto 
entre sí (de ahí el carácter “virtual”, no real, de los trabajos). 


e Indistintamente el SISTEMA EQUILIBRADO puede ser “virtual” y el SISTEMA COMPATIBLE “real”, a 
la inversa, ambos “reales” o incluso ambos “virtuales”. 


e No se presupone el carácter elástico lineal del material, pudiendo existir términos de deformaciones 
de origen térmico o reológico, deformaciones anelásticas, etc. 


Así, el PRINCIPIO DE LOS TRABAJOS VIRTUALES es aplicable para los materiales usuales de edificación en todas sus 
fases de comportamiento mecánico: en el hormigón armado, que exhibe deformaciones anelásticas diversas (fisuración, 
fluencia lenta, etc.); y en el dominio plástico de las piezas metálicas. Para nada se recurre en su demostración a la 
hipótesis de existencia de una función de “energía elástica”, que como función de estado no puede de hecho definirse 
para los materiales no conservativos, sin que ello invalide la aplicación del PRINCIPIO DE LOS TRABAJOS VIRTUALES. 


e Los movimientos deben ser infinitesimales (hipótesis de “pequeños” desplazamientos). 


Seguidamente se enuncian y demuestran, dentro del formalismo del método de equilibrio (variables 
f, s, u, e; matrices H, B), los “teoremas directo y recíprocos” de los trabajos virtuales, de los que 
resultan las propiedades anteriormente expuestas. Asimismo se demuestra la simetría de la matriz de 
rigidez dentro del análisis elástico lineal. 
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Teorema “directo” de los trabajos virtuales: Si se cumplen las ecuaciones de equilibrio £=H-s, las de compatibilidad 
e=Bu, siendo u y e infinitesimales, y H=BY, se verifica la igualdad de los trabajos virtuales externo (We+) e interno 
(Wins). 

Demostración: 
Wes =f!-u=(H-s)'u=(B*.s)u=s!.(B?).u=st.B-u=st.e=W n:, ya que para deformaciones infinitesimales W;¡n¿=s!-e 

Primer teorema “recíproco” de los trabajos virtuales: Si se cumple Wez¿=W:¿n+ para cualesquiera deformaciones 
y movimientos infinitesimales compatibles (e=Bu), entonces las solicitaciones y las acciones verifican f=B*-s 
Demostración: 
fu=Wer=Win+=s'«e=s!-B-u para todo u implica que £t=s?.B, de donde f=B*.s. 

Corolario: H=B*, como se deduce identificando entre sí las expresiones f=B*.s y f=H:s con f y s arbitrarios. 


Segundo teorema “recíproco” de los trabajos virtuales: Si se cumple Wes: =W:;¡n+ para cualesquiera solicitaciones 
y acciones equilibradas (f=H-s), entonces las deformaciones y movimientos, supuestos infinitesimales, verifican 
e=H'u 
Demostración: 
s*.e=Wint=Wex="f'-u=(H-s)'.u=st-H!.-u para todo s implica que e=H!'.u 
Corolario: B=H?, como se deduce identificando las expresiones e=B-u y e=H!-u con u y e arbitrarios. 

Para movimientos no “pequeños”, se cumplen las ecuaciones de compatibilidad e=H*-u=B-u si u y e represen- 
tan movimientos y deformaciones “pequeñas” (infinitesimales) a partir de un estado previo de deformación de la 
estructura. 


Los tres teoremas anteriores se expresan de un modo más compacto como: 
B=H! > Went = Wins 


Simetría de la matriz de rigidez dentro del análisis elástico lineal: Alser D=D*, transponiendo la matriz de rigidez 
K=B*DB se encuentra ella misma: (B*-D-B)'=B*-D*.(B*)'=B'DB 


2.7  PÓRTICOS DESPLAZABLES. 


Se definen como pórticos desplazables aquellos cuyo análisis por el método de equilibrio no puede 
formularse despreciando los desplazamientos (u, v) de sus nudos. Es decir, por definición, los grados de 
libertad de dichos pórticos deben incluir tanto giros como desplazamientos de todos o parte de sus nudos. 


En principio todo pórtico de edificación debe considerarse como “desplazable”: la utilización de un 
modelo de análisis intraslacional para casos que no reunan las correspondientes condiciones de aplicación, 
constituye no sólo un grave error de concepto, sino también la causa de resultados que pueden llegar a 
ser absolutamente erróneos. 


Ya que abordar directamente el caso general de análisis “con desplazamientos” podría dificultar se- 
riamente la comprensión de los métodos de equilibrio, por razones didácticas, se comienza clásicamente 
su estudio por los modelos intraslacionales. 


También por razones didácticas, el reconocimiento de los tipos estructurales a los que son aplicables 
dichos modelos se orienta de manera inductiva, estudiando casos particulares típicos en la edificación 
donde sí son aplicables. Así, se avanza gradualmente en el conocimiento de los tipos específicos más 
importantes de estructuras flectadas. 


Al comenzar el estudio de las estructuras intraslacionales (sección 2.1) se remitió al presente apartado 
para la formulación de un criterio general de aplicación de dichos los modelos. El eventual incumplimiento 
del criterio obliga, lógicamente, a considerar la estructura como “desplazable” a efectos de su análisis por 
el método de equilibrio. 


Este criterio general de intraslacionalidad se puede enunciar de la manera siguiente: 


e Condición “suficiente” de intraslacionalidad: Un pórtico “es” intraslacional si, prescindiendo de los 
voladizos y suponiendo todos los nudos articulados, no se obtiene un mecanismo ni una estructura 
crítica, y si, además, pueden considerarse nulos o despreciables los alargamientos axiles de su barras. 
Normalmente en los pórticos adecuadamente arriostrados puede admitirse el cumplimiento de este 
criterio, salvo la discusión del efecto de la elongabilidad de los arriostramientos que se ilustra 
mediante el ejemplo del apartado 2.9. 
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e Condición “necesaria y suficiente” de intraslacionalidad: Un pórtico puede analizarse como in- 
traslacional si y sólo si, además de ser nulos o despreciables los alargamientos axiles de su barras, 
el equilibrio de fuerzas de los nudos es posible hipotéticamente tomando esfuerzos cortantes nu- 
los en todas sus barras (es decir, haciendo intervenir sólo los esfuerzos axiles N y las reacciones 
exteriores que puedan existir en cada nudo). Si existen cargas aplicadas en el interior de las bar- 
ras (además de posibles cargas en los nudos), para verificar este criterio deben reemplazarse por 
cargas estáticamente equivalentes aplicadas en los nudos extremos de cada barra. En el caso de 
pórticos no arriostrados, o arriostrados de una forma incompleta, debe observarse que ninguna 
de las circunstancias siguientes son, por sí solas o conjuntamente, “suficientes” para asegurar la 
intraslacionalidad: 


— Simetría de geometría, sustentación, constitución y acciones. 
— Inextensibilidad de las barras. 


— Cargas solamente gravitatorias. 


Para demostrar tal afirmación, considérese por un momento su contraria, y tómense como con- 
traejemplos uno o varios de los ejemplos considerados en este texto, donde se analizan pórticos 
que son desplazables pero simétricos (pórtico a dos aguas del apartado 2.10), con barras supuestas 
inextensibles (los pórticos del presente apartado, por ejemplo), y/o con cargas verticales (el mismo 
pórtico a dos aguas antedicho, o el primero de los dos casos del presente apartado). 


Efectivamente, el primer pórtico desplazable que se analiza en este apartado (figura 2.9) se considera 
con una hipótesis de carga “vertical”, a pesar de lo cual se determina un desplazamiento horizontal 
no nulo de sus nudos superiores. 


La hipótesis de inextensibilidad, salvo en el caso de tirantes u otras barras con deformación axil no 
despreciable, se mantiene en general como una hipótesis simplificativa generalmente admisible y muy útil 
para reducir el número de grados de libertad de los modelos de análisis. 


Ello es esencial para poder modelizar, con un número de parámetros reducido, las deformadas de 
pórticos simples clásicos en el estudio de la teoría de estructuras. Con un número de grados de libertad 
reducido son viables desarrollos sencillos, cuya formulación “manual” es indispensable para alcanzar los 
objetivos didácticos que se indican luego. 


Posteriormente se explica que dichas técnicas de reducción de grados de libertad raramente se utilizan 
en el software actual de análisis de estructuras, ya que la capacidad y velocidad de los ordenadores actuales 
permiten cálculos con un elevadísimo número de variables, y resultan en general más sistemáticos y simples 
de programar los modelos que incluyen la totalidad de los grados de libertad. 


Sin embargo se subraya que la utilización de los programas de análisis a modo de “caja negra” no 
puede suplir la etapa de aprendizaje, que de una manera u otra debe pasar por “entrar dentro de la 
estructura”, para asimilar cómo trabajan los tipos estructurales básicos, y qué relaciones de causa a 
efecto se derivan de las diversas condiciones de equilibrio, compatibilidad y constitutivas. 


Sin “saber” estas cuestiones básicas, el manejo del software de análisis de estructuras es tan peligroso 
como conducir un vehículo sin conocer las normas de circulación, como lo demuestra la enormidad de 
los errores de concepto que cada vez proliferan más en los proyectos con cálculos mecanizados de sus 
estructuras. 


Sin embargo sí es cierto que “el valor añadido” de la enseñanza de la teoría de estructuras no incluye 
los aspectos puramente algebraicos o aritméticos del proceso de análisis, por lo que se propone y fomenta, 
su resolución mediante aplicaciones informáticas simples (operaciones de álgebra matricial disponibles en 
pequeñas calculadoras, programas de matemática simbólica, etc.), como forma de “liberar” el análisis 
estructural propiamente dicho de tales aspectos. 


Así, el énfasis de estas lecciones se pone en: 


e Saber plantear todas las condiciones del análisis estructural de cada problema (equilibrio, compat- 
ibilidad y constitutivas). 
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e Saber interpretar críticamente los resultados del análisis, con los diagramas de esfuerzos y defor- 
madas como eficaz forma de expresión de los mismos. 


La formulación adoptada del método de equilibrio se ha elegido precisamente con objeto de separar 
de una manera clara las cuestiones antedichas (sin omitir ninguna de las mismas), de las puramente 
algebraicas O numéricas. 


Pasando a considerar con detalle el caso analizado en la FIGURA 2.9, se observa que la hipótesis de in- 
extensibilidad permite, siempre dentro de la teoría de primer orden (movimientos pequeños, teóricamente 
infinitesimales), suponer que: 


e Los desplazamientos verticales de los nudos a y b pueden omitirse: en efecto, si se desprecian los 
acortamientos de los ejes de los soportes, puede demostrarse que dichos desplazamientos verticales 
son infinitésimos de segundo orden respecto de los desplazamientos horizontales u, y Uu5 . 


e u, YU¿=Ua», puesto que el alargamiento AL,p=u; - U, Se supone prácticamente nulo. 


Así pues, de los seis movimientos posibles: U,, Va, 0, (del nudo a), y us, v», O, (del nudo b), sólo 
se consideran tres incógnitas independientes entre sí: 0.,, Op, Uay ; es decir, se trata de un problema 
analizable con sólo tres grados de libertad efectivos, que constituyen los términos del vector u'=|0,, 
0y, Uas]. Con dichos parámetros queda definida la deformada del estado [ u], que se ha representado 
en la figura con los sentidos considerados como positivos de cada variable (deformada “teórica”). Para 
un desarrollo metódico en materia de signos, es conveniente, como en cualquier problema de análisis, 
plantearlo suponiendo “positivas” las variables (ya sean cinemáticas o estáticas), dejando para el final 
del proceso la verificación de que los signos de los resultados finalmente obtenidos coinciden con los que 
por intuición del comportamiento de la estructura pueden y deben preverse desde el principio. 


Así, en este caso, parece obvio que por efecto de la carga aplicada en la viga, que va a provocar giros 
“hacia el centro” de la misma, los giros 6, y Oy sean respectivamente horario y antihorario, y que, en 
caso de desplazarse los nudos horizontalmente (como de hecho sucede), lo hagan hacia el lado que opone 
menor rigidez (hacia el lado cuyo soporte sea más flexible). Aunque efectivamente se verifica a posteriori 
lo antedicho, no se hace uso de ello a priori en la formulación de las ecuaciones del problema, pues ésta 
es la manera ordenada de plantearlas sin incurrir en errores en los signos de sus términos (lo cual es 
absolutamente esencial, como es lógico). 


Una vez formuladas las variables cinemáticas u*=[9,,, 01, ua»), el criterio de “correspondencia” (aparta- 
do 2.6), según el cual los términos f; y u; se ordenan de manera “correspondiente” desde el punto de 
vista del trabajo mecánico, indica que los términos de £ son m,, m», fa», siendo m,=-qL?/12, m=q1? /12 
(momentos nodales obtenidos igual que en una estructura intraslacional), mientras que fa.» debe obtenerse 
sumando todas las fuerzas “que trabajan” por efecto del desplazamiento u,¿=u¿=U4J), de la planta ab, es 
decir, todas las fuerzas horizontales que pudieran estar aplicadas en los nudos a y b y/o en la barra ab, 
más las posibles reacciones horizontales de los soportes en el estado [f], cuando sobre ellos actúan cargas 
perpendiculares a su directriz (este supuesto se considera en el ejemplo que sigue). 


En este caso f.¿=0; pero resulta conveniente continuar el desarrollo del ejemplo tomando tanto para 
f.y como para ma y my valores sin especificar, reemplazándolos por los datos conocidos f'=[-qL?/12, 
qL?/12, 0] en el momento de resolver las ecuaciones. De esta manera se asegura una vez más no trastocar 
signos de los términos de las mismas, lo que en este caso afectaría a los signos de los resultados u*=[9,, 
0», ua] (como se deduce de la conservación del trabajo virtual W «¿=f*-u), dificultando su interpretación 
final. 


Puesto que no se consideran los alargamientos de las barras, en la expresión del trabajo virtual 
Win:=s!'-e no intervienen los esfuerzos axiles N (cuyo trabajo virtual es Y N-AL=0 al tomarse AL=0 
para todas las barras), por lo que las solicitaciones (s) y deformaciones (e) de las barras son las mismas 
que en los modelos intraslacionales; es decir, en este caso, s*=[my, ..., mg] y e*=[p1, ..., Pe], siempre 
ordenadas de manera “correspondiente” desde el punto de vista del trabajo mecánico. 


La diferencia fundamental con las estructuras intraslacionales es que los giros “relativos” (y) de las 
barras no coinciden en general con los giros “absolutos” (0) de los nudos a los que se unen, como es obvio 
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al observar las deformadas representadas en la figura. La relación de todos ellos se obtendrá luego, al 
escribir las ecuaciones de compatibilidad. 


Por ejemplo, considerando el soporte izquierdo, por giros “relativos” (yp1, (pz) se entienden los que se 
refieren al eje de la barra en su posición final (eje “corrotado”), mientras que los giros “absolutos” de 
los nudos (9=0 en el empotramiento y 0=0, en el nudo a) se refieren a los ejes generales fijos (x, z) de 
la estructura. Por lo tanto las ecuaciones de compatibilidad e=B-u cambian radicalmente con respecto 
a las de una estructura intraslacional. En este caso, al existir un desplazamiento (u,») perpendicular al 
eje de la barra, no coinciden los giros relativos con los absolutos, como puede apreciarse en la figura. 
Naturalmente en los dibujos los movimientos están muy amplificados, pero sus verdaderos valores se 
suponen “pequeños”, lo que justifica el aproximar los ángulos (en radianes) por sus tangentes (así, por 
ejemplo, p¡=Uua,/3.5m). 


El uso de los giros “relativos” (y) permite expresar las relaciones entre momentos y giros de las barras, 
s=D-e, exactamente igual que en el caso de los pórticos intraslacionales, tal como se indica en la figura. 


En efecto, si se toma como eje de referencia de los giros del soporte izquierdo su propio eje en la 
posición final, en el caso de que dichos giros fueran nulos (p1=p2=0), la deformada de la barra sería 
recta (recuérdese que en el estado [u] no se consideran cargas aplicadas en el interior de las barras), y, 
por no existir curvaturas, también serían nulos los flectores M=EI.C. En suma, los momentos extremos 
serían m;=-M(0)=0 y m2=M(3.5)=0 para la barra en cuestión; lo cual es coherente con el resultado de 
aplicar s=D-e, cuando p¡=p2=0. Esto no sucedería si como ejes de referencia se tomasen los ejes fijos 
de la estructura, ya que en el estado citado de deformada recta los giros serían (pi =Uap /3.5m=(2 y por 
lo tanto al aplicar s=D-e se obtendría mi, 40m», resultado absurdo en ausencia de flexión de la barra. 


Para escribir las ecuaciones constitutivas s=D-e manejando en e los giros absolutos, es preciso incluir explícitamente 
los desplazamientos; y, por lo tanto, de acuerdo con Win¿=s!-e, en s aparecen también las fuerzas de extremo de barra. 
Aun siendo esto correcto (de hecho se trata del formalismo más divulgado del llamado cálculo “matricial” de estucturas, 
que puede seguirse en la bibliografía dada en el apéndice C), las dimensiones de las matrices resultan entonces excesivas 
para desarrollos manuales, desvirtuando en carácter sencillo e intuitivo del anterior planteamiento. Por otra parte la matriz 
D resulta entonces singular, denunciando ello que se están mezclando en las ecuaciones constitutivas movimientos de sólido 
rígido además de auténticas “deformaciones”. Pero la condición de equilibrio de cada una de las barras que deben verificar 
las solicitaciones “equilibradas” s permiten descontar en s!-e el trabajo debido a los movimientos de sólido rígido de cada 
barra, con lo cual se encuentra de nuevo que dichos trabajos pueden formularse tomando sólo en consideración los giros 


respecto de los ejes “corrotados”. 


Por las razones expuestas en el apartado 2.6, en las presentes lecciones se considera como objetivo 
docente prioritario saber formular explícitamente las condiciones de equilibrio de cada problema, como 
forma didáctica de continuar, incluso para las estructuras hiperestáticas, trabajando sobre las condi- 
ciones de equilibrio en tanto que herramientas fundamentales del análisis de estructuras, cuya estabilidad 
depende ante todo de que no se vulnere ninguna de dichas condiciones. 


La lectura de los términos de f'=[m,,, m», fas] orienta claramente acerca de qué ecuaciones de equilibrio 
deben escribirse, y en qué orden. Resulta pues un procedimiento muy ordenado y sistemático; además de 
la forma indicada se obtiene una pauta del modo de proceder al enfrentarse por primera vez al estudio 
de un nuevo tipo estructural. 


Se trata de escribir las ecuaciones de equilibrio correspondientes en la forma f=H-s, y extraer seguida- 
mente de ellas la matriz de coeficientes HI como pieza principal de este método de análisis. 


Habiendo estudiado previamente los pórticos intraslacionales, las ecuaciones de equilibrio de momentos 
de los nudos a y b se escriben sin ninguna dificultad, tal como se indica en la figura. 


La única novedad en las ecuaciones de equilibrio radica en la tercera ecuación (correspondiente a fa); 
es decir la que se obtiene expresando el equilibrio de fuerzas horizontales de la planta del pórtico, o “grupo 
cotraslacional” ab que incluye a todos los nudos que participan del mismo desplazamiento U¿=U¿=Uab- 
Según se indica en la figura, resulta una sencilla ecuación de equilibrio Va + V5=f,», donde intervienen 
los esfuerzos cortantes de los soportes, que como se ve son las solicitaciones que equilibrarían la posible 
fuerza horizontal exterior f,»; al particularizar posteriormente f,y=0 deberá resultar obviamente V2=- 
Vs, pero de momento no se hace uso de esta condición (que sin embargo debe verificarse al final, como 
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parte de las condiciones de equilibrio del problema, además de las habituales condiciones de equilibrio de 
momentos de cada nudo). 


Puesto que en el vector s se han tomado como variables independientes los momentos de los extremos 
de las barras, es preciso relacionar con ellos los cortantes Va y Vs. Para ello, el equilibrio de momen- 
tos de cada soporte, permite expresar su cortante V en función de sus momentos inferior y superior, 
reemplazándolos luego en la ecuación de equilibrio de fuerzas horizontales. 


Se obtienen tres ecuaciones en las que intervienen como incógnitas sólo los momentos my, ..., mg de 
los extremos de las barras. La matriz de coeficientes, de 3 filas y 6 columnas, es la matriz de equilibrio 
H, indicada en la figura. 


La relación general B=H!* demostrada en el apartado 2.6 permite obtener automáticamente la matriz 
de compatibilidad B por una simple transposición de filas por columnas. 


Por lo tanto no es necesario duplicar el trabajo de análisis deduciendo ahora explícitamente las ecua- 
ciones de compatibilidad e=B-u. Pero sí es importante saber interpretarlas, es decir, saber escribir las 
diversas filas de e=B-u (ver figura) y verificar su coherencia geométrica sobre el dibujo de la deformada 
de la estructura. De esta manera: 


e Se comprueba, que la matriz B, y por lo tanto también la H, son correctas. 


e 5e visualiza cómo se mueve y deforma el pórtico objeto de este ejemplo. 


Lógicamente es perfectamente válido partir de las ecuaciones de compatibilidad e=B-u, extraer de 
las mismas la matriz B, y obtener H por transposición (o bien, directamente, formar la matriz de rigidez 
K=B'DB), siempre que se compruebe, dada su absoluta importancia, que las expresiones f=B*.s son 
las ecuaciones de equilibrio correctas del problema. 


Puesto que, por lo indicado, la matriz constitutiva D es la misma que la conocida: de los modelos 
intraslacionales, y sigue siendo válida la cadena de sustituciones f=B*-s=B*.D-e=B*.D-B-u (EQU-ECO- 
COM), resulta el sistema de ecuaciones típico del método de equilibrio, K-u=f, de donde se obtienen 
los movimientos u=K”7*-f, y a partir de ellos, los esfuerzos s=D-B-u del estado fu], que se completan 
sumándoles algebraicamente los esfuerzos sf del estado [f], de manera análoga a lo ya puesto en práctica 
para estructuras intraslacionales. 


Nótese que en todo este proceso, después de “almacenados en memoria” el vector de acciones nodales 
f y las matrices H y D, se opera exactamente igual que para estucturas intraslacionales. Para resolver el 
problema, basta manejar como “objetos” sucesivamente las matrices: 


B=H';, K=H-D-B; u=K”! +; s=sp; 4+D-B -u 


de las cuales, salvo a efectos de comprobación (como por ejemplo la matriz B, según lo ya indicado), basta 
solamente obtener explicitamente u y s, manejándose el resto de matrices simbólicamente (módulos op- 
eraciones de álgebra matricial disponibles en pequeñas calculadoras, programas de matemática simbólica, 
etc.). 


De esta manera se trabajan explícitamente los conceptos básicos del análisis estructural (equilibrio, 
compatibilidad, ecuaciones constitutivas) y se libera de aspectos tediosos y que no forman parte propia- 
mente del mismo, como son las operaciones puramente numéricas, sin caer en el manejo de “cajas negras” 
informáticas. 


Una vez que esta técnica resulta familiar, los tiempos de análisis son muy cortos (y se centran en 
lo realmente importante), de manera que puede y debe dedicarse un tiempo, fundamentalísimo, para la, 
comprobación de los resultados. Como mínimo, se deben realizar siempre las sencillas comprobaciones 
siguientes: 


e Equilibrio de momentos de cada nudo 
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Equilibrio de fuerzas horizontales de cada planta (para pórticos ortogonales de varias plantas; planta 
única en el ejemplo analizado). 


Equilibrio global del pórtico. 


Equilibrio de fuerzas de cada barra en direción de su eje y en dirección perpendicular. 


Carácter equilibrado del diagrama de flectores de cada barra 


Propiedades gráficas de los diagramas de esfuerzos, estudiadas en la teoría de vigas. 


Verificación del carácter “lógico” de los movimientos calculados, según consideraciones similares a 
las que para este ejemplo en particular se expusieron al principio de este apartado. 


Coherencia de las curvaturas de las barras (C) en la deformada final, con los signos de los flectores 
(M) en cada punto, según se desprende de la relación C=M/El. 


En la programación en ordenador de esta técnica de análisis, análogamente a lo indicado en apartados 
anteriores, las matrices B y D se dividen en cajas correspondientes a cada una de las barras, obviando 
áreas extensas de dichas matrices que para grandes estucturas están constituidas por ceros. La matriz 
de rigidez de la estructura (KK) se obtiene efectuando para cada barra los productos del tipo BD¿B, y 
sumando o “ensamblando” entre sí las matrices así obtenidas. 


Puede fácilmente relacionarse el formalismo aquí expuesto del método de equilibrio con otros más 
clásicos que pueden seguirse en la bibliografía mencionada en el apéndice C (métodos de Cross, “slope- 
deflection”, “cálculo matricial”, etc.). La lectura de esa bibliografía se recomienda enfáticamente como 
medio para adquirir la cultura propia de un universitario, en el sentido de comprender la enorme y rica 
diversidad de aproximaciones posibles a los distintos problemas, y aprender a relativizar enseñanzas que 
de otro modo prodrían parecer dogmáticas. 


El método propuesto, que se sigue desarrollando en los apartados siguientes, se ha elegido por su 
operatividad, su economía conceptual (el mismo formalismo es útil para multitud de problemas, incluso de 
posibles futuros cursos de especialidad), pero sobre todo porque se manejan ineludiblemente los conceptos 
básicos de equilibrio y deformación de las estructuras analizadas, lo cual se pierde totalmente si los cálculos 
se reducen a la utilización más o menos eficiente de una herramienta informática “cerrada”. 


En la figura 2.10 se desarrolla un segundo ejemplo de análisis de un pórtico desplazable ortogonal de 
una sola altura, que permite generalizar algunos aspectos del análisis efectuado en el caso anterior. 


En ambos ejemplos, análogamente a lo ya realizado en el estudio de las estructuras intraslacionales, 
cuando existen cargas aplicadas en puntos intemedios de las barras, el principio de superposición, válido 
dentro del análisis lineal, permite descomponer el estado real de la estructura en otros dos estados: 


e Un primer estado de nudos empotrados [f], donde se aplican las cargas de las barras. 


e Un segundo estado de nudos libres [u], con acciones sólo en los nudos 


A diferencia de las estructuras intraslacionales, en el primer estado es preciso considerar tanto las 
fuerzas como los momentos de empotramiento, y en este segundo ejemplo ello resulta relevante para 
obtener la acción nodal f., equivalente sobre la planta ab. 


Para ello, simplemente tanto las fuerzas como los momentos de empotramiento de las barras se aplican 
cambiados de signo a los nudos de la estructura, y se suman algebraicamente a las fuerzas y/o momentos 
directamente aplicados en los mismos, para obtener las acciones nodales equivalentes de la siguiente etapa. 


Otro aspecto que difiere del primer ejemplo consiste en la existencia de un extremo de barra articulado, 
cuyo giro (ys) puede no considerarse como una variable independiente (puesto que, como se demostró en el 
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apartado 2.2, está relacionado con el otro giro de la barra en la forma pg=-p5/2). El giro pg no se incluye 
en el vector e, que en este segundo ejemplo pasa a tener sólo 5 términos. El momento mg=0 tampoco se 
incluye en el vector s. Análogamente a lo que ocurría al considerar las barras con un extremo articulado 
que formaban parte de estructuras intraslacionales, debe adoptarse la rigidez equivalente 3El/L, en lugar 
de la matriz “4-2-2-4”. 


Según lo ya conocido, como alternativa se podría haber tomado para la barra articulada del segundo 
soporte la matriz general “4-2-2-4”, incluyendo tanto pg como mg en e y s respectivamente; pero entonces 
también el giro 0 del apoyo debería añadirse en u como un grado de libertad más, y por lo tanto sería 
necesario aumentar a 4 el número de ecuaciones de equilibrio, incorporando la correspondiente ecuación 
de equilibrio de momentos del nudo. 


Finalmente, ambos ejemplos ponen de manifiesto que en pórticos desplazables no se cumple la sencilla 
regla de la “transmisión” de momentos a los empotramientos de los soportes en la cimentación (basta, 
para cerciorarse de ello, comparar los momentos obtenidos en dichos extremos con los que resultarían 
de aplicar dicha regla). En pórticos desplazables, dichos momentos deben obtenerse como un término 
más del vector s=sf + D-B-u, tal como puede seguirse en los desarrollos de ambos ejemplos que se 
incluyen en las figuras. Este detalle es sumamente importante en los pórticos de edificación, porque 
típicamente los valores absolutos máximos de los momentos de los soportes se producen precisamente en 
sus empotramientos en cimentación. 


2.8 PÓRTICOS ORTOGONALES DE VARIAS PLANTAS. 


Se generaliza el formalismo del método de equilibrio expuesto en la sección anterior, en la que sólo 
se consideraron estructuras modelizables con un único desplazamiento además de los giros de los nudos, 
a estructuras con un número mayor de desplazamientos, como son las estructuras de edificación más 
comunes: los pórticos ortogonales de varias plantas. 


El objetivo prioritario de estas lecciones, formular explícitamente los distintos sistemas de ecuaciones 
(equilibrio, compatibilidad y relaciones elásticas entre esfuerzos y deformaciones) necesarios en la resolu- 
ción de cualquier problema hiperestático, limita la aplicación manual del método a estructuras de unas 
pocas plantas como la que sirve de ejemplo en este apartado (FIGURA 2.11). 


El estudio de estructuras ortogonales porticadas con un elevado número de plantas y/o vanos se deja así 
para una lección posterior a las contempladas en el texto actual. Las adecuadas reducciones de grados de 
libertad, por condiciones de repetitividad, simetría y/o antimetría, permiten su análisis mediante modelos 
con pocas variables: modelos para los que puede utilizarse manualmente el formalismo propuesto en estas 
lecciones, preservándose el objetivo docente indicado. Naturalmente la mecanización en ordenador del 
método de equilibrio amplía enormemente el número de grados de libertad de las estructuras analizables. 


Puesto que sigue aceptándose la hipótesis de inextensibilidad de las barras (para una discusión detalla- 
da de esta hipótesis véase el apartado 2.9) todos los desplazamientos de los nudos de las estructuras que se 
estudian en este apartado son horizontales; sin embargo el desarrollo que se expone sería totalmente válido 
si además de estos fuese necesario considerar las componentes verticales de los movimientos. Bastaría 
añadir los términos adecuados al vector f, de tal modo que el producto u'-f evaluase convenientemente 
Wez: (según lo expuesto en la sección 2.6), y en consecuencia habría que considerar las ecuaciones de 
equilibrio de fuerzas correspondiente a los nuevos elementos del vector de acciones nodales; puede seguirse 
un ejemplo al respecto en el apartado 2.10. 


También se utiliza por primera vez el concepto de ANTIMETRÍA como alternativa para reducir signi- 
ficativamente el número de grados de libertad de algunos modelos estructurales. Se hará referencia a los 
resultados de la sección 2.2 relativos a la matriz equivalente D.gu de una viga con deformada antimétrica. 


Examinando la deformada de la FIGURA 2.11, es claro que las vigas de este ejemplo responden al caso 
particular p,=3 considerado en el apartado 2.2: los giros de los extremos izquierdo y derecho de cada 
viga coinciden en valor y en signo. 
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Para cualquier barra que cumpla esta condición, se tiene: 


_ ABI BL, _8El 
mi = L Pi E Pa = L P1, 


como ya se mostró en el apartado indicado. Esta relación es idéntica a la que existiría si se situase úna 
articulación (ARTICULACIÓN VIRTUAL) en el punto medio de la viga. Para una barra de longitud L/2 con 
su extremo articulado, la rigidez equivalente según lo expuesto en el apartado 2.2 es precisamente: 


3EI  6El 
HZ" E 
Por este motivo en el ejemplo de la FIGURA 2.11 se puede analizar sólo media estructura. 


El razonamiento anterior se basa en la noción de REPETITIVIDAD de las variables cinemáticas de la estructura (en este 
caso, la de los giros izquierdo y derecho de cada viga). 


La simplificación del cálculo que ello permite (en caso de poder aplicarse) es tan notable, que a veces se aplica cuando se 
pretende obtener resultados de carácter sólo aproximado (no exactos). Aunque no se cumplan rigurosamente las condiciones 
p1i=p2 en todas las barras se introducen estas articulaciones virtuales. Esta es la técnica utilizada para desarrollar los 
MODELOS APROXIMADOS para el análisis de pórticos típicos de edificación mediante SUBESTRUCTURAS. 


Otra forma, basada en consideraciones físicas, de llegar al modelo indicado, se apoya en la ANTIMETRÍA 
(anti-simetría, o simetría negativa) del ejemplo. Observando la deformada de la estructura se aprecia 
que el punto central de las vigas no experimenta ningún movimiento vertical y su curvatura es nula 
(puesto que los movimientos y la curvatura de los situados a izquierda y derecha del mismo tienen signos 
contrarios). Recordando la relación entre flectores y curvaturas se deduce que el flector en dicho punto 
debe de ser nulo. Así, el modelo puede simplificarse sustituyendo las vigas por unas barras ficticias de 
longitud L/2 y con un apoyo deslizante “virtual” en el eje de antimetría del pórtico. 


A la igual conclusión se llega considerando el flector como una pareja de momentos (dorsal y frontal) en las dos secciones 
de cada rebanada de las vigas. Cuando ésta coincide con el eje, sí M0, se tendrían dos momentos (el dorsal y el frontal) 
simétricos, en contra de la hipótesis de antisimetría. Por reducción al absurdo se deduce que debe ser M=0 en el eje medio. 
El mismo razonamiento permite deducir que N=0 en el eje medio. En cambio, suponer un cortante V0 no es contradictorio 
con la hipótesis de antimetría, ya que la fuerza cortante dorsal es antisimétrica de la frontal (“una sube y la otra baja”). 


No debe extraerse la consecuencia de que el análisis de pórticos ortogonales desplazables puede efec- 
tuarse siempre suponiendo antimetría (del mismo modo que las diversas simplificaciones por simetría 
sólo pueden utilizarse en las condiciones adecuadas y con sumo cuidado para no incurrir en gravísimos 
errores conceptuales y numéricos). La simplificación por antimetría sólo es válida cuando se cumplen las 
condiciones siguientes: 


e Simetría geométrica, de sustentación y de constitución. 


e Acciones antimétricas, en el sentido luego indicado, como lo son las parejas de fuerzas horizontales 
P a cada lado del eje de simetría que pueden sustituir a la fuerzas 2P de la figura. 


La antimetría es pues una propiedad específica de algunas estructuras como la de la FIGURA 2.11 (igual 
que ocurría con la simetría), pero es cierto que es frecuentemente “exportable” a muchos problemas de 
pórticos sometidos a acciones horizontales. Se ha elegido el ejemplo como “representativo”. 


La antimetría se define como el producto o encadenamiento de una transformación de simetría y una 
inversión de signo entre dos semi-estructuras: a toda magnitud (fuerzas, momentos, flechas, giros, etc.) 
de una mitad de la estructura, le corresponde en la otra mitad NO SU MAGNITUD SIMÉTRICA, SINO SU 
CONTRARIA. 


La elección de los elementos de f es inmediata en virtud del criterio de “correspondencia”, una vez 
dibujada la deformada de la estructura y definidas las variables cinemáticas (los giros 0, y 0, junto a 
los movimientos horizontales de los dinteles u, y us), que se agrupan, como siempre, en el vector u. En 
la FIGURA 2.11 se han tomado como términos del vector de acciones nodales los momentos m, y my de 
los dos nudos del modelo antimétrico y las fuerzas horizontales f, y fp que actúan sobre cada uno de los 
dinteles. De este modo el producto f!-s representa el trabajo W e: - 
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No se ha representado el estado [f] por ser trivial para el ejemplo propuesto. Al no existir acciones 
sobre las barras de la estructura todas las reacciones, solicitaciones, deformaciones, etc. de dicho estado 
son nulas. Las componentes del vector f son, por lo tanto, directamente las acciones aplicadas sobre 
los nudos, siguiendo las pautas generales expuestas en las secciones previas. Así, m,=mp=0, pero por 
las razones también ya expuestas (garantizar que no se cambian inadvertidamente los signos de alguna 
ecuación y consecuentemente de algún resultado) no se hará uso de estas condiciones hasta el momento 
de resolver el sistema de ecuaciones final. Del mismo modo el desarrollo se realiza tomando dos acciones 
horizontales genéricas fa y fo. 


Debe notarse que en el modelo antimétrico adoptado se sustituyen las acciones “reales” aplicadas en 
los nudos situados a la izquierda de la estructura por unas acciones antimétricas, para lo cual se “traslada” 
la mitad de cada una de ellas al nudo frontal de cada dintel. Esta argucia, lícita para hacer uso en esta 
etapa del análisis de la simplificación por antimetría, debe ser reconsiderada en la fase final del proceso 
en la que, como se ha insistido, debe verificarse el equilibrio de cada nudo con las acciones “reales” que 
sobre él actuan. Con esta salvedad, el vector de acciones nodales que se adopta es £'=[0, 0, P, P]. 


Las ecuaciones que expresan las condiciones de equilibrio de momentos de los nudos a y b son similares 
a las de cualesquiera de los ejemplos hasta aquí tratados. Asimismo la ecuación que expresa la condición 
de equilibrio de las fuerzas horizontales de la planta superior es idéntica a la expuesta en el apartado 
anterior, salvo que sólo interviene en ella un cortante, pues como es patente en la FIGURA 2.11 en el 
modelo antimétrico sólo se considera un soporte en cada planta. 


Obsérvese que la articulación virtual de cada viga se comporta a todos los efectos como un apoyo 
deslizante, en el cual la reacción sólo puede tener componente vertical, por lo que en ningún caso inter- 
vendrá a la hora de plantear la condición de equilibrio según el eje horizontal. La citada componente 
vertical de la reacción de la articulación virtual no será nula, pues en otro caso el cortante en las dos vigas 
cancelaría, lo que a su vez implicaría flector constante y por tanto nulo: único modo de que se cumpla 
a la vez la condición de antimetría. Obviamente la reacción -acción- de media estructura es la acción 
-reacción- de la otra media, y representa por lo tanto el esfuerzo cortante que existe en dicha sección de 
la estructura “real”. 


La única novedad en las ecuaciones de equilibrio se encuentra al plantear la de la planta intermedia. 
En ella intervienen los cortantes de los soportes situados por encima y debajo de la planta, no sólo los 
inferiores como ocurría al expresar la condición de equilibrio de la planta superior (o las del único dintel 
de las FIGURAS 2.9 y 2.10). 


Nuevamente hay que extremar las precauciones en materia de signos, siendo vital la distinción entre 
esfuerzos y fuerzas; de igual modo que al plantear las ecuaciones de equilibrio de momentos se insistió en 
la diferencia entre flectores y momentos. Las incógnitas consideradas son los cortantes de los soportes, 
que como se ha reiterado es conveniente suponer positivas, mientras que la ecuación que se escribe es 
de equilibrio de fuerzas. Así se justifica que las dos incógnitas Va y V3 aparezcan en la ecuación de 
la figura con signos contrarios. Al plantear la condición de equilibrio según el eje horizontal, mientras 
debe tomarse la parte superior (frontal) de la primera rebanada del soporte superior, ha de considerarse 
la inferior (dorsal) de la última rebanada del soporte inferior. En definitiva f-Va+-V3=0, y tras dejar 
en el lado izquierdo de la igualdad el término corespondiente al vector de acciónes nodales (de manera 
coherente con el sistema ecuaciones de equilibrio escrito de la forma f=H-s) resulta f =V32-V3. 


Las dos ecuaciones que expresan la condición de equilibrio de las plantas se han escrito como función 
de los cortantes de los soportes. Puesto que las componentes en el vector s son los momentos m, en los 
extremos de las barras, deben sustituirse todos V; por las expresiones que los relacionan con los momentos 
en los extremos de los soportes. El proceso puede seguirse en la figura y es similar al expuesto en ejemplos 
anteriores. 


Una vez que las incógnitas de las cuatro ecuaciones de equilibrio son únicamente las seis componentes 
del vector s, la matriz del sistema (de cuatro filas y seis columnas) es la matriz H de equilibrio. 


La matriz de compatibilidad B se obtiene por simple trasposición de la anterior según se dedujo en 
el apartado 2.6. 


Aunque parezca superfluo escribir ahora las ecuaciones de compatibilidad, interpretarlas sirve para 
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verificar las ecuaciones de equilibrio y profundizar en el conocimiento de los movimientos y las deforma- 
ciones de las estructuras. 


Las ecuaciones que expresan las condiciones de compatibilidad de los giros de las vigas (filas quinta 
y sexta de la matriz B de la figura) son tan triviales como las de las estructuras intraslacionales, por lo 
que no requieren nuevos comentarios. Las que expresan las condiciones de compatibilidad de los giros 
relativos del soporte inferior (filas primera y segunda de e=B-u) son idénticas a las de los pórticos de una 
sola planta estudiados en el apartado 2.7. La única novedad se encuentra en las ecuaciones que expresan 
la compatibilidad de giros de los extremos del soporte de la planta superior, que indican que los giros 
relativos son la suma del giro del nudo correspondiente (0, o 0) más el giro del eje del soporte. Este 
se estima, en primera aproximación, como el cociente entre movimiento relativo de la planta superior 
respecto de la inferior y la altura del soporte (tg(a = «: en radianes para ángulos pequeños), es decir 
(u¿-Ua)/Ha. 


El razonamiento anterior sugiere tomar como variable cinemática independiente no el movimiento absoluto de la planta 
superior sino el relativo entre ambas plantas. De este modo las ecuaciones de compatibilidad (y por ende de equilibrio) se 
simplifican notablemente. Este camino se explora unas líneas más abajo. 


Ya se ha justificado que las matrices constitutivas de las vigas son las D.¿ correspondiente a una barra 
con deformada antimétrica deducidas en la sección 2.2, mientras que para los soportes deben considerarse 
la matriz estándar de 2x2 según se indicó en la sección anterior. 


A partir de este punto, conocidas las matrices H, B y D y el vector de acciones nodales, el proceso 
es idéntico al seguido para cualesquiera de los tipos estructurales ya estudiados: 


K=B'.D-B, f=K-u, u=K"?-f, s=sf+D-B-u 


Obtenidos los momentos en los extremos de las barras pueden dibujarse los diagramas de flectores, a 
partir de estos los de cortantes y finalmente, por equilibrio de fuerzas de los nudos los de axiles. 


Finalmente han de realizarse las sencillas comprobaciones que se señalan en la seccion 2.7 (equilibrio 
de momentos en los nudos, equilibrio global, etc.), prestando especial atención en este caso al equilibrio 
de fuerzas de cada nudo. Como se indicó al inicio del apartado para poder hacer uso de las condiciones 
de antimetría se traslado la mitad de las acciones de los nudos a y b a sus simétricos. Sin embargo la 
condición de equilibrio de los nudos debe, obviamente, verificarse con las acciones que actuan “realmente” 
sobre ellos. Así se entiende el axil en cada uno de los dinteles, que en rigor para un modelo antimétrico 
debía ser nulo según se justificó previamente (nótese que el valor de este axil es precisamente -P, lo que 
explica que sea lícito trasladar la mitad de la acción de los nudos a y b a sus simétricos, ya que hipótesis 
de inextensibilidad garantiza que la presencia o ausencia de dicho axil no modifique Wins). 


Ya se ha mostrado en casos anteriores cómo modelizar una misma estructura de distintas formas: es 
obvio que en las estructuras simétricas o antimétricas podría no hacerse uso de estas condiciones, que 
en las estructuras con barras con un extremo articulado puede adoptarse un grado de libertad adicional 
asociado al giro del extremo en cuestión en lugar de utilizar la apropiada rigidez equivalente, o que, en 
definitiva, podrían asignarse tres grados de libertad (el giro y en general dos movimientos) a cualquier 
punto interior de la estructura. La discusión que sigue muestra un modo más sutil de obtener modelos 
diversos de una estructura. En lugar de modificar las variables cinemáticas que se consideran, se varía la 
medida que de ellas se realiza, o en otras palabras, se cambia el sistema de referencia respecto del que se 
miden. 


Siguiendo la observación realizada al interpretar las ecuaciones de compatiblidad de la estructura de 
la FIGURA 2.11, en la FIGURA 2.12 se analiza la misma cambiando la variable cinemática us, que medía 
el movimiento del dintel superior referido al sistema global de referencia, por u',, que mide el mismo 
movimiento pero tomando como sistema de referencia uno que se traslada con el dintel inferior. Es decir, 
u', representa el movimiento relativo del dintel superior respecto del inferior. 
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La primera fila de la FIGURA 2.12 muestra las diferencias entre ambas descripciones cinemáticas para la 
estructura analizada. Se indican diversos modos de deformación de ambos modelos. Nótese que mientras 
en el modelo ya estudiado el estado definido por el vector de movimientos u*=[0, O, 1, 0] representa una 
deformada en la que el dintel superior permanece en la situación original, un vector de movimientos con 
iguales componentes para el nuevo modelo u*=[0, 0, 1, 0], representa una configuración radicalmente 
distinta de la estructura. En este caso ambas plantas tiene un movimiento “global” común. El estado 
u'=[0, 0, 1, 0] del primer modelo corresponde al u'*=[0,0,1,-1] del segundo y el estado u'*=[0,0,1,0] del 
nuevo modelo con la descripción original se escribiría como u*=[0,0,1,1]. 


Esta nueva interpretación cinemática de la estructura debe afectar al vector de acciones nodales. Como 
se probó en la sección 2.6 el producto u'-f, o el u'*-f', debe evaluar el trabajo de las fuerzas exteriores, 
que obviamente será el mismo para ambos modelos. 


Como es habitual el criterio de “correspondencia” permite identificar inmediatamente los términos 
del nuevo vector f'. Nada nuevo hay que decir respecto a las componentes ma, mp y fo, que serán las 
mismas del primer modelo. Respecto a la componente f',, se obtiene sin dificultad recordando que debe 
contener todas las fuerzas que trabajan cuando se produce el movimiento correspondiente del vector u!. 
Es claro que u*=[0,0,1,0] moviliza todas las acciones horizontales situadas sobre el dintel inferior, siendo 
por lo tanto f',=f, +f+. 


A igual conclusión se llega en la figura mediante un ejercicio meramente algebraico. Se reescribe el 
Wez: del primer modelo en función de las variables cinemáticas del segundo, sin más que introducir el 
cambio de variables u¿=u', y U4=u'¿+u'», tal como se indica en la figura. 


Sirva lo anterior de ejemplo clarificador de la afirmación del apartado 2.6 según la cual la expresión 
Wes = -u = )) f;:u; permite generalizar el concepto de “acción” f; a cualquier magnitud estática, o 
conjunto de ellas, que “trabajan por” un cierto movimiento u;. 


El cambio en el vector de acciones nodales conlleva necesariamente el de las ecuaciones de equilibrio 
que se deben plantear. Así, para el nuevo modelo en lugar de la ecuación de equilibrio de fuerzas del 
dintel de la planta baja aislado, debe escribirse la del grupo cotraslacional formado por toda la estructura 
situada sobre dicho dintel. De este modo se da una interpretación desde la perspectiva estática (que se 
prima en estas lecciones) del cambio de variables cinemático que originó el nuevo modelo. 


La formulación final se realiza en la figura sobre la base de las ecuaciones de equilibrio, obteniéndose 
del modo habitual la matriz H. La matriz de compatibilidad B se obtiene por transposición de H y la 
interpretación de sus filas es inmediata a la luz de los párrafos anteriores. 


Como se previó al discutir las condiciones de compatibilidad del primer modelo las matrices de com- 
patibilidad y de equilibrio son más sencillas, pero a cambio los términos del vector f son más elaborados. 
De hecho, la necesidad de incluir en f', las acciones del propio dintel y las de la planta superior requiere 
un momento de reflexión, frente al carácter más sistemático del desarrollo del primer modelo. 


Como se ve los vectores s y e no se modifican, por lo que tampoco lo hace la matriz D. 


El resto del proceso es idéntico al ya descrito varias veces. La única diferencia se encuentra al 
inspeccionar el nuevo vector de movimientos. No podía ser de otro modo, pues su cuarta componente 
representa ahora el movimiento relativo de la planta superior, como se verifica en la figura. 


La relación entre los dos modelos pueden interpretarse, como se insinúa en la exposición, como un 
simple cambio de variables. Así, mediante el procedimiento general para cualquier cambio de variables 
que se enunciará en el apartado 2.10 se puede reproducir de un modo elegante y conciso la discusión 
anterior. 


Se concluye este apartado en el que, entre otros aspectos, se trata el problema de la antimetría, recor- 
dando que cuando el eje de simetría de la estructura coincide con un soporte la simplificación propuesta 
deja de ser válida (de igual forma que en los problemas simétricos el planteamiento es radicalmente dis- 
tinto si sobre el eje de simetría existe o no un soporte). En este caso la reducción por antimetría requiere 
sustitutir el soporte por uno de rigidez mitad. 
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2.9 PÓRTICOS ARRIOSTRADOS. 


Los objetivos del presente apartado son: 


e Introducir en el formalismo desarrollado en las lecciones anteriores, donde se han considerado es- 
tructuras formadas exclusivamente por piezas flectadas, barras biarticuladas solicitadas sólo por 
esfuerzos axiles (N). Las modificaciones en el proceso de análisis seguido en apartados anteriores 
son muy sencillas y permiten ampliar los tipos estructurales analizables a una amplia variedad: 


l 


Pórticos ortogonales de una o varias plantas, arriostrados en todas o parte de ellas. 


Estructuras atirantadas o colgadas mediante tirantes extensibles. 


Estructuras flectadas en general en las que no se puedan considerar nulos o despreciables los 
alargamientos axiles de su barras. Este caso se explica luego, en el presente apartado. 


| 


Pórticos no ortogonales, en los que, independientemente de si admiten o no la hipótesis de 
inextensibilidad, ésta puede llegar a constituir, no una simplificación del modelo, sino más bien 
todo lo contrario. Para un estudio de esta cuestión, puede seguirse el pórtico no ortogonal 
analizado en el apartado 2.10 


e Discutir la eficacia de los arriostramientos. Normalmente en los pórticos adecuadamente arriostrados 


puede admitirse que la condición de intraslacionalidad del apartado 2.5, en el cual no obstante se 
remitía a la presente lección para el estudio del efecto de la elongabilidad de los arriostramientos. 


Es obvio el interés práctico de esta cuestión, ya que la caracterización de un pórtico como in- 
traslacional permite utilizar los modelos de análisis mucho más sencillos, ya estudiados en lecciones 
anteriores. 


También esta discusión es fundamental para cuestiones que trascienden del ámbito puramente lineal 
de este conjunto de lecciones, pero que se incluyen en el mismo curso de Análisis de Estructuras de 
Edificación: pandeo de soportes; métodos de cálculo plástico. A dichos efectos puede “exportarse” 
la caracterización o no de un pórtico arriostrado como intraslacional a que se llega en el presente 
apartado, ya que, a pesar de efectuarse en contexto elástico lineal, es perfectamente suficiente para 
dilucidar si en un pórtico dominan la rigidez de flexión de sus piezas o la rigidez de conjunto de sus 
arriostramientos, ante acciones horizontales diversas. 


El ejemplo desarrollado en la figura 2.13 es muy simple, pero suficiente para los objetivos indicados. 


Obviamente, se trata del mismo pórtico analizado en la figura 2.10, al que se le han incorporado 


dos diagonales; de este modo se pueden comparar los esfuerzos y movimientos ante las mismas acciones 
horizontales, y evaluar la mayor o menor incidencia de las diagonales añadidas. 


Las diagonales se suponen de muy reducida rigidez a flexión frente al valor EI/L indicado en la figura, 


que corresponde al resto de las barras del pórtico. Por este motivo, sólo se considerará eficaz la diagonal 
traccionada, pues la diagonal comprimida pandear para un esfuerzo |N| muy reducido que se supone no 
altera el balance de fuerzas del modelo. 


Este criterio (considerar sólo la familia traccionada de diagonales de arriostramiento, según el sentido 


supuesto de las acciones horizontales, que normalmente pueden actuar en los dos sentidos posibles), es 
habitual en el cálculo de arriostramientos, pero requiere un comentario de sus condiciones de aplicación: 


e El pandeo de la diagonal comprimida debe producirse en el dominio elástico, y por lo tanto, sin 
deformaciones plásticas irreversibles 


e El pandeo debe realmente producirse para compresiones muy pequeñas, de tal forma que al pandear 
la barra no puedan inducirse fenómenos dinámicos alarmantes. 
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Puede suponerse que las condiciones anteriores se verifican cuando la esbeltez (A) de las diagonales es lo 
suficientemente elevada, aunque sin sobrepasar lo límites normativos. Naturalmente, además, la diagonal 
comprimida no debe ser estrictamente “necesaria” para equilibrar las acciones exteriores (ya que en tal 
caso se produciría el colapso de la estructura); pero ello se cumple merced a existencia de la otra diagonal, 
traccionada, gracias a lo cual el pandeo de la diagonal comprimida se desarrolla por “acortamiento 
impuesto” limitado por el resto de la estructura. Así dicha diagonal comprimida simplemente abandona, 
su directriz rectilínea de forma casi imperceptible, o por lo menos de manera reversible. 


Comparando los desarrollos de las dos figuras 2.10 y 2.13, se observa que las variables nodales f'=[m,,, 
mo, fas] y u'=[0., 0,, ua] son idénticas; la expresión del trabajo virtual “externo” W.,¿=f!-u no se ve 
pues formalmente alterada. 


Sin embargo, en el vector de solicitaciones (s) es preciso incluir la tracción de la diagonal, Na (donde 
el subíndice 4 hace referencia a dicha barra, la cuarta de la estructura). 


En efecto, Na obviamente va a intervenir al formular la ecuación de equilibrio de fuerzas horizontales 
de la planta, que constituye la tercera fila del sistema f=H-s escrito en la figura. En la misma, el esfuerzo 
Na de la diagonal se incluye como último término de s. 


En la ecuación de equilibrio citada, tal como se indica en la figura, aparece la componente horizontal 
de Na, lo cual determina la aparición del coeficiente cos(a.) en la posición correspondiente de la matriz de 
equilibrio H. Con respecto a la matriz de equilibrio del apartado 2.7, se ha añadido una séptima columna. 
cuyos términos son 0, 0 y cos(a). 


En el apartado 2.2 ha quedado clara la importancia del principio de los trabajos virtuales en la 
fundamentación del método de equilibrio, por lo que para preservar la expresión del trabajo virtual 
interno Wins=s'-e, cualquier término que se incorpore en s debe corresponderse con su término dual en 
e, y viceversa. 


Del estudio de las estructuras trianguladas que precede (entre otras) a estas lecciones, se toma la 
expresión del trabajo virtual interno correspondiente a una barra articulada cualquiera: N¿-ALgz, siendo 
AL; el alargamiento de la barra en el estado compatible. En este caso tanto N¿ como ALy4 corresponden 
al mismo estado “real” de la estructura, aunque como sabemos ello no es necesario en general en el 
formalismo de los trabajos virtuales. 


Así pues, en este ejemplo se tiene: 
s'=[m,, ..., m6, Na) 
e'=[p1, ..., Ps, ALa] 
Wini=s!-e=m1 -p1+...+m6po+N4-ALa 
Al transponer la matriz de equilibrio, B=H*, la última columna añadida en H se corresponde con 


una nueva fila (la última) en la matriz de compatibilidad, tal como consta en la figura. Si se denomina 
Ba a dicha fila correspondiente a la diagonal (cuarta barra de la estructura); se tiene: 


B4=[0, 0, cos(a)] 
Por lo tanto, tomando la última de las ecuaciones de compatibilidad e=B-u, se puede escribir: 
AL¿=B4-u=0-0,+0-0,+cos(a)-Ua=ucos(a) 
Esta ecuación de compatibilidad admite la sencilla interpretación geométrica (ver figura) de que el 


alargamiento AL4 debe coincidir con la proyección sobre el eje de la diagonal del desplazamiento del nudo 
b al cual se une, ya que su otro extremo está fijo (apoyo inferior izquierdo). 


Ello sirve, entre otras cosas, como comprobación de que son correctas la fila y columna añadidas 
respectivamente en B y H, con respecto a las matrices análogas del apartado 2.7. También permite 
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verificar una vez más la armonía entre las condiciones de equilibrio y compatibilidad en todo formalismo 
“bien temperado” de análisis estructural. 


En las ecuaciones entre esfuerzos y deformaciones s=D-e, a las relaciones entre momentos (m) y giros 
relativos (p) de las barras consideradas en el apartado 2.7, es preciso añadir aquí la relación entre -los 
nuevos términos de s ye, es decir: 

_ EA, 


N, 
4 La 


+ ALa, 


como se deduce de la expresión conocida AL¿=N¿ L¿/EA;z del alargamiento de una barra biarticulada 
cualquiera. 


Por lo tanto, en la matriz constitutiva D es preciso añadir la constante EA4/La en la posición corre- 
spondiente de su diagonal principal, tal como se indica en la figura. 


Esto completa el conjunto de modificaciones del algoritmo de análisis con respecto al formulado en el 
apartado 2.7. El resto del proceso es formalmente idéntico: 


K=B'.D-B, u=K”! £, s=s+D-B-u 


También es posible ensamblar la matriz de rigidez de la estructura mediante la expresión ya conocida 
K=>K;, siendo K¿=B;*-D;-B; las contribuciones de cada barra, tal como se indica en la figura. 


Para cualquier estructura flectada en la que se considere necesario o conveniente considerar los 
alargamientos axiles de todas o parte su barras, se puede proceder análogamente, incluyendo los co- 
rrespondientes esfuerzos N, en el vector s, los alargamientos AL; en el vector e, y las constantes EA 4 /Lz 
en las posiciones correspondientes de la diagonal principal de la matriz constitutiva D, como en el ejemplo 
aquí expuesto, aunque los esfuerzos Ny pertenezcan a barras flectadas cuyos momentos mz también estén 
incluidos en s. 


En efecto, la expresión AL;=Nz L¡,/EAz del alargamiento de una barra biarticulada, es aplicable 
también a cualquier barra (flectada o no) solicitada por un esfuerzo normal N¿ constante, como se 
deduce mediante la cadena de igualdades AL=e7-L=(00/E)L=(N/A)-L/E=NL/EA, donde por claridad 
se ha omitido el subíndice k de la barra, y se ha adoptado la expresión del alargamiento unitario del eje 
baricéntrico de la barra, eo =N/EA, válida dentro del análisis elástico lineal. 


Un ejemplo de este proceso puede seguirse en el estudio del pórtico no ortogonal del apartado 2.10 


Como límite, con el presente modelo de análisis es resoluble cualesquier estructura articulada, isostática 
o hiperestática. 


Si todas las barras están biarticuladas, el vector de esfuerzos incluye solamente los esfuerzos normales 
Ny, el vector de deformaciones solamente contiene los alargamientos AL¿, y la matriz constitutiva D 
tiene como únicos coeficientes no nulos las rigideces axiles EA¿/Lz, en su diagonal principal. 


Para una estructura articulada, la matriz H se deduce a partir del sistema de ecuaciones de equilibrio 
de fuerzas f. y f, de cada nudo, y aunque su formación es muy fácilmente programable, en este curso 
básico de análisis se pide saber escribir dichas ecuaciones de equilibrio una por una (tanto en la fase de 
planteamiento inicial del mismo como en la comprobación final de los resultados, que como mínimo debe 
incluir de nuevo la verificación de todas las condiciones de equilibrio del problema). 


Volviendo al ejemplo del pórtico arriostrado de la figura 2.13, en la misma se ha incluido un parámetro 
adimensional 


que se utiliza para relacionar la rigidez de los dos sistemas rigidizantes ante acciones horizontales “macla- 
dos” en esta estructura: solicitaciones axiles, especialmente la del tirante; y solicitaciones de flexión del 
resto de las barras. 


Para simplificar en primera aproximación dicho estudio, se considera el caso de que la viga sea infini- 
tamente rígida. Entonces 0,=0=0», y de la tercera fila de f=K-u se deduce: 
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£,=K31:04+K32:05+K33 Us =K33 Uas 
Por lo tanto, para el caso de ser la viga infinitamente rígida, el coeficiente K33=fap/Uay mide la rigidez 
ante acciones horizontales. 


Desarrollando la contribución de la barra diagonal a la matriz de rigidez de la estructura: 
K,=B4*-D¿-Ba 


(donde, según lo ya expuesto anteriormente, B4=[0, 0, cos(a)] y D4 =[EA4/La)), se encuentra la matriz 
K, expresada en la figura, cuyo único término no nulo es: 


el cual se ensambla precisamente en la posición (3,3) de K 
Procediendo análogamente con el resto de las barras, se encuentra su contribución a K3g : 


12- El 12- El 
dildo: a a E 

donde la sumatoria se extiende únicamente a los soportes, y representa la rigidez de flexión (es decir, la 
rigidez del pórtico no arriostrado) para el caso de una viga infinitamente rígida. La rigidez del pórtico no 
arriostrado real será lógicamente inferior, dependiendo de la mayor o menor rigidez de la viga. Así, en el 
ejemplo del apartado 2.7 se obtuvo fas /Uay=80 KN /0.01392 m=4310 kN /m, mientras que 24-El/L3=9600 
kN/m (algo más del doble), rigidez que por las razones luego indicadas resulta conservadora a los efectos 
de evaluar la eficacia del arriostramiento. 


Continuando con el caso de una viga infinitamente rígida (0,=0=04) se tiene: 
fa=Kpr UK pr Up =for Hr, 


donde for=KppUap es la parte de fa, absorbida por el arriostramiento, y fpr=K pr-Ua» la parte de f.¿ que 
produce flexiones en el pórtico. Despejando Uy = y reemplazando en ff, =K nu), Se encuentra: 


como es lógico, ya que la hipótesis simplificativa 0,¿=0=0, reduce el número de grado de libertad a uno 
(u,»), en cuyo caso, tal como se vio en la sección 2.3 (estructuras intraslacionales de un grado libertad) 
y se acaba de volver comprobar, la acción exterior se reparte proporcionalmente a las rigideces de los 
diferentes elementos o sistemas involucrados. 


Como la expresión aproximada K,,=)>312-El/L3 sobrevalora K ;,, resulta efectivamente conservadora 
a efectos de estimar la parte de las acciones horizontales no absorbida por el arriostramiento (fp). 
Reemplazando los valores anteriores de K y, y Kyr, se puede pues escribir: 


E, A 
fu < TE ESP 
siendo B=(EA4/La)/(El/L*) en el ejemplo, y en general: 


as EAarr/Larr 
N EII > 


donde para pórticos regulares de varias plantas, la sumatoria se extiende a todos los soportes de una 
planta típica o representativa, y EA arr/Larr corresponde también al arriostramiento de la misma planta. 


B 


Un criterio clásico en construcción metálica, recogido en los Eurocódigos, y obviamente razonable, 
consiste en considerar un pórtico como arriostrado si el arriostramiento reduce sus desplazamientos hor- 
izontales en, al menos, un 80 % . Según esto, el pórtico de la figura 2.13 NO podría considerarse como 
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adecuadamente arriostrado, ya que el desplazamiento obtenido es de 3.27 mm, que supera a 0.2-13.92 
mm=2.78 mm, siendo 13.92 mm el desplazamiento del pórtico sin arriostrar (apartado 2.7). 


Sin necesidad de analizar completamente los dos pórticos (el arriostrado y el no arriostrado), una 
forma conservadora de verificar dicho criterio consiste en limitar: 


de $ e a 0-2, 

ab 1+ 24 

de donde B:cos*(a) >96, o bien: B-cos?(a) >100, CONDICIÓN CONSERVADORA PARA PODER CONSID- 
ERAR “ARRIOSTRADOS” LOS PÓRTICOS REGULARES de una o varias plantas (pero, en principio, no 
más de 10). Como comprobación puede verificarse que tampoco se cumple en el ejemplo analizado 
(B-cos*(a)=33<100). 


2.10  PÓRTICOS NO ORTOGONALES. 


Para exponer cómo puede aplicarse el método de análisis estudiado en estas lecciones a los pórticos 
no ortogonales, se considera el ejemplo de la figura 2.14, correspondiente a un pórtico a dos aguas, 
biarticulado. Aunque se trata de un tipo particular de pórtico no ortogonal, muestra la manera apropiada 
de incluir en el análisis piezas inclinadas que pueden aparecer en otros tipos de pórticos muy diversos. 


Además del objetivo indicado, el ejemplo que se analiza permite las siguientes enseñanzas de diseño 
estructural, “exportables” a otros casos que pueden plantearse: 


e La directriz de la pieza no adopta, obviamente, la forma antifunicular de las cargas (que sería en este 
caso una parábola de segundo grado), pero hasta cierto punto puede considerarse una aproximación 
burda a dicha forma, haciendo concesiones a aspectos constructivos tales como el de conformar las 
fachadas y la cubierta: con superficies planas, y constituir la estructura a base de piezas rectas. 
En este sentido, el pórtico a dos aguas puede considerarse como un “arco imperfecto”, en el que 
aparecerán flexiones, por no coincidir su directriz con la línea de presiones. (En realidad eso sucede 
en mayor o menor grado en todo arco real, ya que la línea de presiones puede siempre variar 
por alternancias de cargas variables, y porque el eje de la pieza no es fijo, sino que experimenta 
desplazamientos que pueden producir flexiones de segundo orden a veces muy importantes, como 
en todo sistema comprimido). De todas formas, si se compara el diagrama final de flectores de este 
problema con el que correspondería a un pórtico simple ortogonal, se aprecia una notable reducción 
de los flectores máximos, si el pórtico está suficientemente peraltado. Relacionar el problema 
con conceptos de estática gráfica que se suponen previamente conocidos, permite interpretar y 
comprobar el diagrama final de flectores del dintel del pórtico, cuya forma debe coincidir con la 
diferencia de ordenadas de un arco parabólico (con arranques en los apoyos articulados) y la directriz 
poligonal realmente adoptada. 


e La razón de ser de las articulaciones en cimentación se comprende bien si se analiza el mismo 
pórtico con empotramientos en su base, lo que se deja como ejercicio. El cociente e=M/N en las 
secciones inferiores de los soportes (que es cero en el pórtico biarticulado), adopta entonces valores 
considerables si las luces son medias-altas y las cargas de cubierta son ligeras, que es precisamente 
el campo típico de aplicación este tipo estructural. Puesto que el cociente “e” es la excentricidad 
de la reacción sobre la cimentación, se encuentra que el posible beneficio que se obtendría para 
la estructura, al hacerla más coaccionada en su base, se anula por el inconveniente de requerir 
cimentaciones desproporcionadamente anchas por razones de su estabilidad al vuelco. En cambio 
en el pórtico biarticulado las únicas reacciones son N (compresión) y V (cortante en la base del 
soporte) que son perfectamente asumibles si las reacciones horizontales se pueden anclar en la 
solera, o en su defecto si se equilibran mediante un elemento traccionado que una las bases de los 
dos soportes (adecuadamente protegido contra la corrosión, por quedar enterrado). Nótese que por 
equilibrio general de la estructura las dos reacciones horizontales son autoequilibradas (de igual 
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EA 
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b | m2 0. Ma q 
Empotramiento : mi Ua Toa P2 ma 
deslizante 5 v, f > 3 m3 
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valor absoluto y sentidos contrarios); únicamente las acciones horizontales originan una resultante 
horizontal neta, generalmente de menor entidad que las anteriores reacciones individuales de cada 
soporte debidas a las cargas verticales. 


e Una tercera articulación (por ejemplo, en la clave), convertiría la estructura en isostática, lo. que 
no impediría la aplicación de la misma técnica de análisis aquí expuesta, pero ofrecería alternativas 
más simples al quedar todos los esfuerzos y reacciones estáticamente determinados. Por compara- 
ción con los diagramas de flectores de los pórticos triarticulados, que se suponen estudiados con 
anterioridad a estas lecciones, puede observarse que los valores máximos son bastante similares. 
Pero en cambio el hiperestatismo reduce los movimientos de la estructura (lo cual es en general una 
constante del diseño de estructuras de edificación), y este aspecto puede ser determinante para el 
dimensionado de este tipo de estructuras de luces medias-altas, ya que las flechas y desplazamientos 
en general aumentan con la luz mucho más rápidamente que los esfuerzos de flexión. Además, con la 
rótula en clave se obtiene en general más una complicación que una simplificación constructiva, por 
lo que el pórtico triarticulado se emplea raramente, como en general los pórticos isostáticos, salvo 
circunstancias especiales (por ejemplo, asientos excepcionalmente elevados que no puedan combat- 
irse con cimentaciones profundas, ya que los asientos no provocan solicitaciones si la estructura es 
isostática). 


La discusión anterior deja claro que el pórtico biarticulado es una estructura de grado de hiperestatismo 
igual a uno, por lo que es un tipo estructural apropiado para analizarlo por el método de compatibilidad, 
lo cual se deja como ejercicio, como alternativa al método de equilibrio o “de la rigidez” aquí expuesto. 


El análisis de este pórtico por el método de equilibrio se detalla en la figura 2.14, donde se siguen las 
diferentes etapas, muy sistemáticas, características de este método ya estudiado en las lecciones anteriores. 
Unicamente deben recordarse o destacarse algunos aspectos específicos, lo que se efectúa a continuación. 


En la etapa [f], al estar las cargas referidas no a la longitud L de las piezas inclinadas, sino a su 
proyección horizontal L. (tal como suele hacerse de hecho en la práctica) los esfuerzos de empotramiento 
perfecto coinciden con los de una hipotética viga empotrada horizontal de longitud L.,, como ya se vio 
en el ejemplo de la figura 2.7 


El modelo adoptado en la figura hace uso.de la simetría del problema para suponer que el desplaza- 
miento horizontal de la clave u, es cero, así como su giro 0,, pero en modo alguno el desplazamiento 
vertical v,, cuya omisión constituiría un gravísimo error de concepto y trastocaría dramáticamente to- 
dos los resultados. Tomar v,=0 sólo estaría justificado si existiera un apoyo o un soporte (supuesto 
inextensible) que impidiera dicho desplazamiento vertical (Véase la discusión final del apartado 2.5). 


Se plantea pues en la figura 2.14 un modelo de tres grados de libertad u*=[0,,, ua, vo]. Posteriormente 
se discute acerca de la conveniencia o no de aplicar la hipótesis de inextensibilidad del dintel para reducir 
el número de grados de libertad. Esta discusión se establece con carácter más académico que práctico, ya 
que la posible aplicación o no de la hipótesis de inextensibilidad (que, como toda hipótesis simplificativa 
es optativa, no obligatoria), en el caso de piezas inclinadas generalmente complica más que simplifica el 
análisis, por lo que su adopción carece en general de sentido práctico. 


Sin embargo la supuesta inextensibilidad del soporte sí se ha aplicado, ya que el desplazamiento 
vertical v, no se ha incluido en el modelo. En las aplicaciones informáticas, generalmente suele ser más 
sistemático incluir los movimientos u, v, 0 de todos los nudos. En este ejemplo, no obstante, considerar 
va € AH en el modelo (que podrían, perfectamente, incorporarse) redundaría en una cierta pérdida de 
sencillez de la exposición, la cual se ha considerado prioritaria. 


Como siempre, el número y carácter de las variables cinemáticas orienta acerca del número y naturaleza 
de las ecuaciones de equilibrio que deben formularse, y que se corresponden también con las acciones 
nodales f'=[ma, fra, £,4], que en la figura se han calculado siguiendo la regla general de aplicar al nudo 
los momentos y fuerzas de empotramiento de la(s) barra(s) cambiados de signo, sin olvidar añadir las 
posibles acciones directamante aplicadas a los nudos que pudieran existir. 


Como se ve al escribir las ecuaciones de equilibrio, en las mismas aparece el esfuerzo Nay del dintel, 
que por lo tanto debe incluirse en el vector de solicitaciones s*=[m,, m2, m3, Na»], puesto que éste se 
obtiene precisamente reuniendo todos los esfuerzos de las barras que intervienen en aquéllas ecuaciones. 
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La matriz de equilibrio H se obtiene simplemente formando una matriz de tantas filas como ecuaciones 
de equilibrio y tantas columnas como términos de s, y asignado a los términos de la matriz los coeficientes 
de cada término de s en cada una de las ecuaciones citadas, de manera que éstas adopten la forma general 
f=H:s 


El vector de deformaciones es e*=[g1, 2, P3, AL], como se deduce de sumar los trabajos m; -p1+ma* 
(2 -+m3 - (3 +Nas - AL e identificar la expresión obtenida con Winy=s!-e. 


Las ecuaciones de compatibilidad son e=B-u, obteniéndose B por transposición de H. En la figura se 
interpretan geométricamente estas ecuaciones, que posteriormente se estudian también con más detalle. 


Para acabar de resolver el problema, simplemente basta con escribir la matriz constitutiva D, tomando 
para el soporte la rigidez 3El,/H (barra con extremo articulado), para los momentos ma y mg del dintel 
la matriz estándar de flexión 4-2-2-4, y para Nay la rigidez axil EA,/L, tal como se dedujo en el apartado 
anterior. 


Las operaciones matriciales habituales: 
K=B'.D.-B=K?*f s=sp+D-B-u 


seguidas de la representación de diagramas de esfuerzos y de la deformada, completan el proceso de 
resolución, que se detalla en la figura. Como siempre, deben realizarse comprobaciones de coherencia de 
la deformada y de equilibrio general de la estructura y de sus nudos y elementos, algunas de las cuales se 
indican en la misma figura. 


Los resultados obtenidos confirman las previsiones inicialmente formuladas en este apartado, tanto 
en lo relativo a la forma del diagrama de flectores como en que respecta a los elevados valores de los 
desplazamientos, que respetan justamente los límites menos exigentes aplicables en estructuras de edifi- 
cación (estos límites se dan para estructuras metálicas en la norma española EA95 y en el Eurocódigo 
3). Los desplazamientos calculados serían totalmente inaceptables en caso de existir elementos construc- 
tivos rígidos dañables por deformación. Aún bajo dicha hipótesis, la especificación en el proyecto de una 
contraflecha de ejecucion apropiada sería en este caso indispensable. 


Este mismo ejemplo sirve para introducir un concepto importante para cualquiera de los métodos 
de análisis que por uno u otro procedimiento desembocan en un sistema de ecuaciones del tipo K-u=f. 
Se trata del posible “mal condicionamiento” de la matriz de rigidez K, Puesto que dicho formato es 
hoy el estándar en el software de análisis de estructuras, puede afirmarse que casi cualquier programa 
de análisis puede fallar para estructuras “mal condicionadas”, y que en determinadas circunstancias el 
programa puede no detectar que se están obteniendo resultados incorrectos. 


Ello constituye un ejemplo de las muchas razones por las cuales una intuición clara del comportamiento 
y forma de trabajo de las estructuras, junto con las diversas comprobaciones cualitativas y cuantitativas 
del análisis, sobre las cuales tanto se ha insitido en estas lecciones, son absolutamente necesarias para 
que el empleo de los programas de análisis de estructuras no pueda entrañar un riesgo similar al del uso 
de un vehículo sin el conocimiento de las más elementales reglas de conducción. 


Este peligro de “mal condicionamiento” es consustancial a la implementación del método de la rigidez 
en ordenador digital, y se produce cuando la matriz de equilibrio es cuasi singular, det(KX) x= 0, de tal 
forma que las rutinas de solución del sistema de ecuaciones K-u=f no produzcan una señal de error, aun 
cuando sus resultados sean incorrectos. 


Esta situación puede encontarse cuando en el ensamblaje de la matriz K se combinan rigideces enorme- 
mente dispares de los diferentes elementos que integran la estructura. En tal caso, como el número de 
digitos significativos que se almacenan en memoria para cada uno de los coeficientes K;¿ es, naturalmente, 
limitado, puede suceder que se pierdan las contribuciones de uno o varios elementos a un mismo coefi- 
ciente K;¿ si son mucho menores en valor absoluto que una de ellas, de tal forma que debieran sumarse a 
cifras decimales que más allá de la posición duodécima o decimosexta, sencillamente no se almacenan. Sin 
llegar a este extremo, también puede suceder que aun formándose correctamente todos los términos de K, 
sin embargo los efectos de uno o varios de los elementos invulocrados se perdieran en las manipilaciones 
artméticas necesarias para que el ordenador resuelva el sistema de ecuaciones K-u=f 


El mal condicionamiento es particularmente probable si un elemento o constante de rigidez muchísmo 
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más elevada que el resto no se une directamente a la cimentación, sino a través de otros elementos 
mucho más flexibles. "Todo usuario de software de análisis de estructuras debería conocer el riesgo de mal 
condicionamiento en tal caso. 


Por ejemplo, en el pórtico de la figura, la rigidez axil EA/L del dintel puede ser numéricamente muy 
superior a las rigideces de flexión del dintel y/o del soporte. Al efectuar el producto B*-D-B, aparecen por 
ejemplo para el soporte términos del tipo 3El,/H* ; dividiendo EA/L por el anterior valor, se encuentra 
el cociente adimensional siguiente: 

A,H y» 
3AJL” ” 


siendo A=H/i,=esbeltez del soporte en el plano del pórtico (i¿=radio de giro=y/I,/As. En el ejemplo 
A? es superior a 1000, lo que muestra que las rigideces axil y de flexión pueden ser muy dispares entre 
sí. Para órdenes de 1000 y aún muy superiores es improbable que se produzca mal condicionamiento de 
la matriz de rigidez del sistema. Pero si se adoptase una rigidez EA/L ficticia muy superior a la real, se 
estaría incumpliendo la regla de “no conectar directamente a la cimentación rigideces muy elevadas” (la 
rigidez axil del dintel) a través de otras mucho más flexibles (la rigidez de flexión del soporte), con un 
riesgo claro de mal condicionamiento numérico. 


Por lo tanto no debe pretenderse en el problema de la figura 2.14 simular la “hipótesis de inextensi- 
bilidad” del dintel introduciendo datos ficticios EA/L anormalmente elevados. 


Respetando esta sencilla e importante regla, el procedimiento de análisis indicado en la figura resul- 
ta más apropiado y sencillo que la introducción de la “hipótesis de inextensibilidad” del dintel. Ésta 
puede derivarse como una ligadura de los desplazamientos ua y Vy impuesta por la cuarta ecuación de 
compatibilidad (cuarta fila de e=B-u) tomando en la misma AL “exactamente” igual a cero: 


AL = 4 : cos(a) — uy - senla) = 0 
Ua = 0) :tg(a) 


Esta condición permite considerar u, como variable dependiente de vs (o al revés), y por lo tanto 
reducir a dos el número de grados de libertad efectivos, en lugar de los tres iniciales. 


Además de lo anterior, el esfuerzo Na» se elimina de las ecuaciones de equilibrio simplemente sumando 
la tercera más la segunda multiplicada por tg(a), como puede comprobarse. 


Para mantener en todo momento la simetría entre la formulación estática y cinemática del problema, 
es conveniente plantear las ligaduras de la forma siguiente: 


e Escribir el cambio de variables cinemáticas: u=L-u”, siendo en este ejemplo: 


ul == [0., La, Ub), w' = [0.,, vs) 


1 0 
L=| 0 tg(a) 
o 1 


como se comprueba simplemente efectuando el producto matricial L-u' 
e Reducir el número de ecuaciones de equilibrio en la forma: 
f=L'.f=L'-H:s 
Con ello, lo que se está haciendo en este caso simplemente es eliminar Na, de las ecuaciones de 


equilibrio de la manera antes indicada sumando la tercera más la segunda multiplicada por tg(a), 
como puede comprobarse directamente realizando el producto matricial L*-H-s. 


*pe9.191] 9P SOPe13 Sp UPTIESUMPUOO) :GT'Z PINE 


0 —2ente) al 


pr 
A 
cosía 
pa A 1-2 sente) E 
3 pe: 
e: 1 4 0 de 


AL 0 cos(a) —sen(a) 
Última fila: 


AL = u, - cos(a) — vs - sen(a) 


Variables 
Nudos Barras 
mi 
0% m + 
$e f A e ma 
( zb Pa m3 


AL=0 


o vo-tg(a) 


0 pa 


Interpretación geométrica 


po 


AL=0=uz-cos(a) — v - sen(a) pa] [p-=0 aa we 
2 ¡ 1 _senla) _ cos(a) Si 
Us = 05 :tg(a) |= e Es y, - tg(a) 
3 | $ + 0 ás 
b 
Cambio de variables AL] lo cos(a) —sen(a) 
as LA 1.0 a e 030 pe : y 
= <= e => = u 
oy=tl jo ted lA P2 = 1 AO] Ml =P". u' 
u,=0; :tg(a) PB P3 1 7 
u= L w LAL] lo 0 B"=B- L 
EQU, ECO COM Solución 
H= wo 4 E o | KX'-B'*.D'.B' u=K'-1.f= 
Ale A. 6) Aa 0.005012 
coslaJL  cos(aJL H B'=¡* Zeosía) 0.06268 
K'=B"* «.D-B” Ñ 
tgla 
SEl 2El q 17 | 
L L L i s'=D'-B'u'= 
n_ |2EL 4ELn y K'=L*-B'.D-B-L 1109 
D'=|Z L É 
SEL Bu=| B ] 80.81 
0. =i o] H1—[26820 1506 120.8 
—|-1506 599.1 


DIAGRAMA DE M (kN m) 


DESPLAZAMIENTOS (m) 


67e-1 


0.6268e- 


Acciones nodales 


Wez=u! -f 
=04 - Ma + Un “fea +0 * f2b 
=0, Ma + -tg(a) Si Fa +0%-* fas 
=0, Ma +05 - (tg(a) * fza + $25) 
=u”- [ma,tg(a) - fza + fal! 


f'=[mo, tg(a) > Íza + fas] 
f=Lt-f 


Cambio de variables alternativo: 
Se toma u, como parámetro independiente 


. 
= 47 A 
Vb = ly» Ya = Ue, 6%=0a4 


021 [1 0 
8, "Et. 

ua|=|0 1 41] qbrE 

1 La 
Ub 0 E = Ma 

ee z Plis 
u-= L al 

Acciones nodales 
Wezi=u! -f 


=02 Ma + Ua fza +0 f2b 
=0, Ma +Uo - fza + za Les 


=0, Mag + Un * (fea =$ des) 
=u*. [ma, Ízo + des] 


£=[ma, feo + tala J 
f=Lt.f 


v9 
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Es inmediato comprobar que f'=L*-f es el cambio “dual” de variables estáticas que mantiene in- 
variante el trabajo virtual W ¿1 =f'-u=f?-u!. 


Al no intervenir ya Nay en las ecuaciones de equilibrio, se puede quitar del vector de solicitaciones, 
con lo cual éste se reduce: s“=[m,, ma, ms] 


Las ecuaciones de equilibrio quedan así: =H'-s', siendo en este ejemplo: 
0 1 1 


1 
n= A 1 ta(a) 
Lcos(a)  Lcos(a) H 


Al ser s*=[m,, m», mg], la matriz constitutiva se reduce a una matriz 3x3 donde intervienen sólo 
las rigideces de flexión: 


o 
1 | 281 4Ér 
71% 7 

H 


e Realizar las operaciones matriciales habituales con las nuevas matrices y variables: 


B'= H" K = p” -D'. B'u' = Ko L A y = s, +D'-B'.u 


Aunque interesante como ejercicio teórico, es obvio que este proceso es menos simple que el planteado 
inicialmente, por lo que resulta poco práctico. Solamente para cálculos puramente manuales, puede ser 
más ventajoso resolver el sistema de dos ecuaciones K'.u'=f" que el sistema inicial de tres ecuaciones 
K.u=f. Lo mismo es generalizable para sistema con muchos más grados de libertad, para los que la 
capacidad y rapidez de los ordenadores actuales ha convertido en obsoletas muchas técnicas de reducción 
de grados de libertad, que retienen un valor más académico que práctico. 


El mayor interés de las expresiones F=L*-f, u=L-u', K'=L*-K-L es que constituyen un procedimiento 
perfectamente general para cualquier cambio de variables en el análisis de estructuras, y de hecho engloban 
como caso particular las ecuaciones estándar del método de la rigidez, f=B*. s, e=B-u, K=B*.D-B, 
que pueden así considerarse como una tranformación entre variables estáticas y cinemáticas y rigideces 
“internas” (s, e, D) y “externas” (f, u, K), con B y B* como “matrices de paso”. 
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MÉTODOS DE EQUILIBRIO 
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VALORES DE INTEGRALES QUE CONTIENEN PRODUCTOS. 


Apéndice B. 
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Figura B.1. Vigas hiperestáticas simples. Momentos de empotramiento. 


( creeme Nm Porra 


AR 


7 
“Mag 


-M2E 


-M2E 


70 


VIGAS HIPERESTÁTICAS SIMPLES. MOMENTOS DE EMPOTRAMIENTO. 


Apéndice C. 


BIBLIOGRAFÍA PARA CONSULTA. 


Análisis de secciones y elementos: 
Gere, J.M. y Timoshenko, S.P.: “Mecánica de Materiales”. G.E.I., 1986 


e Resumen de estática gráfica. Estructuras Trianguladas. Métodos clásicos de cálculo de pórticos: 
Timoshenko, S.P., y Young, J.M. : “Teoría de las Estructuras”. Ed. Urmo, 1975 


Métodos de cálculo elástico: 


Livesley, R. K. : “Métodos matriciales para cálculo de estructuras”. Ed. Blume, 1970. 
West, H.H.: “Analysis of Structure”. Wiley, 1989. 


e Métodos de cálculo plástico: 
Baker, J., Heyman, J.: “Plastic Design of Frames”. Cambridge, 1969 


Colección de ejercicios resueltos: 
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